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Abstract 

Let G be a fixed connected multigraph with no loops. A random n-lift of G is 
obtained by replacing each vertex of G by a set of n vertices (where these sets are 
pairwise disjoint) and replacing each edge by a randomly chosen perfect matching 
between the n-sets corresponding to the endpoints of the edge. Let Xq be the 
number of perfect matchings in a random lift of G. We study the distribution 
of Xq in the limit as n tends to infinity, using the small subgraph conditioning 
method. 

We present several results including an asymptotic formula for the expectation 
of Xq when G is (i-regular, d > 3. The interaction of perfect matchings with short 
cycles in random lifts of regular multigraphs is also analysed. Partial calculations 
are performed for the second moment of Xq, with full details given for two example 
multigraphs, including the complete graph K±. 

To assist in our calculations we provide a theorem for estimating a summation 
over multiple dimensions using Laplace's method. This result is phrased as a 
summation over lattice points, and may prove useful in future applications. 

Keywords: random graphs, random multigraphs, random lift, perfect matchings, Laplace 's method. 
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1 Introduction 



Throughout, let G be a fixed connected multigraph with g vertices and no loops. For 
simplicity we assume that V(G) = [g] := {1, . . . ,g}. A random n-lift of G is a random 
graph on the vertex set V\ U V% U • • • U V g , where each Vi is a set of n vertices and 
these sets are pairwise disjoint, obtained by placing a uniformly chosen random perfect 
matching between Vi and Vj, independently for each edge e = ij of G. Denote the 
resulting random graph by L n {G). The perfect matching corresponding to the edge e of 
G is called the fiber corresponding to e, which we denote by F e . Note that the degree 
of v G Vi in L n {G) is equal to the degree d G {i) of vertex i in G. In particular, if G is 
(i-regular, then so is L n (G). We are interested in asymptotics as n tends to infinity 

This model of sparse random graphs was introduced and studied in a series of papers 
by Amit, Linial, Matousek, and Rozenman (2j [3j HJ [12] . Linial and Rozenman [T2] studied 
the existence of a perfect matching in L n (G) and described a large class of graphs G 
for which L n (G) a.a.s. contains a perfect matching (for n even, at least). This class 
contains all regular graphs and, in turn, is contained in the class of graphs having a 
fractional perfect matching (see Section [3] for a definition). Observe that if G has a 
perfect matching then every lift of G has at least one perfect matching. 

In this paper we study the number of perfect matchings in L n (G) in the limit as 
n tends to infinity, where G is a graph with a fractional perfect matching. To do this 
we use the small subgraph conditioning method, which provides a concentration result 
based on the second moment method conditioned on the number of small cycles. For a 
concise description of the method, see [TH Theorems 9.12 and 9.13]. 

Let Xq be the number of perfect matchings in L n (G). To apply the small subgraph 
conditioning method, asymptotic expressions for EXg and K(X G ) must be found. Then 
the limit of the ratio K(X G ) / (KXg) 2 is compared against a quantity which depends upon 
the interaction of perfect matchings and short cycles in L n (G). 

In Sections [3] and H] we write the first and second moments of Xq as multiple sums 
of some explicit terms, and then estimate the sums by Laplace's method. This is a 
standard method for similar moment estimates, and in particular, it has been used in 
several papers on random regular graphs. (See for example [TTJ, Chapter 9] and the 
references given there.) However, in the present paper, each summation is over an index 
set of rather high dimension with a number of side conditions on the indices, while in 
many previous applications the summations are only over one or two variables. To assist 
with these calculations, we present a general theorem (Theorem 12.31) that encapsulates 
Laplace's method for a general situation, with sums over a lattice in a subspace of M. N . 
We do this both because we think that it clarifies the argument in the present work, 
and because we hope that it might be useful in future applications. The necessary 
terminology and notation is introduced in Section [2j where Theorem 12.31 is stated. The 
proof of Theorem 12.31 can be found in Section 

Using this machinery we prove an asymptotic formula for EXg for any connected 
regular multigraph G with degree at least three (see Theorem 13.61) . However, two diffi- 
culties (one algebraic and one analytic) have prevented us from obtaining an asymptotic 
formula for K(X G ) in the same generality, though we have partial results in Theorem 14.21 



2 



and Lemma 14.31 We illustrate these results by calculating K(Xq) for two multigraphs: 
specifically, for the complete graph K± and for the multigraph consisting of two vertices 
and three parallel edges, which we denote by K\. These calculations were performed 
with the aid of Maple. (A file containing the Maple commands is available from [20J.) 

In Section [5] we prove the necessary results relating to short cycles in random lifts 
( Lemmas 15 . 1 [ [5721 and Corollary 15.41) . As corollaries, using [Tij Theorem 9.12] we obtain 
a concentration result for Xq in our two illustrative examples (see Corollaries 15.71 and 

ESD- 

One of the most interesting questions on random lifts is the problem of existence of 
a Hamilton cycle. There is a conjecture (attributed to Linial) that a random lift of K± is 
a.a.s. hamiltonian. Indeed, we believe that a.a.s. L n (G) is hamiltonian for all connected 
c?-regular loop-free multigraphs G with d > 3. (This is known to be true when G is a 
multigraph with exactly two vertices and at least three edges: see Remark [1.11 below.) 
Burgin, Chebolu, Cooper and Frieze [B] showed that a.a.s. L n [K g ) is hamiltonian when g 
is large enough (see also [7] for the directed case). The arguments in [5] are combinatorial 
and utilize the celebrated idea of Posa. For small g, we feel that the small subgraph 
conditioning method may be a fruitful line of attack, as it has been very successful for 
studying Hamilton cycles in random regular graphs (Robinson and Wormald [UJ, [18] . 
see also [TT1 Chapter 9]). This remains an open problem. 

Remark 1.1. We allow the multigraph G to have multiple edges. The simplest case is 
when G consists of only two vertices, with d parallel edges between them. The random 
lift L n (G) then is a random bipartite (multi)graph obtained by taking the union of d 
independent random matchings between two sets of n vertices each. Such sums have 
been studied in [15], where they were shown to be contiguous to random bipartite d- 
regular (multi)graphs. The latter, in turn, is known to be a.a.s. hamiltonian (see [16] 
for a standard, second moment method proof). Hence for this small multigraph G with 
d > 3, the random lift L n (G) is a.a.s. hamiltonian too. 

Remark 1.2. Random lifts of multigraphs with loops can also be formed. As in [2], the 
fiber corresponding to a loop is given by the n edges ia(i) for a random permutation o 
of [n]. This is a random 2-regular (multi)graph, denoted by P(n) in [TTJ Remark 9.45]. 
While we do not allow loops in our current work, for several reasons, we believe that 
the results here can be extended to multigraphs with loops. A simple and interesting 
case is when G consists of a single vertex with d/2 loops (d even). Then L n (G) consists 
of the sum (union) of d/2 independent copies of P(n). Such sums have been shown to 
be contiguous to random d- regular (multi)graphs in [8]. 

Acknowledgements. This work was partly done at the Isaac Newton Institute, Cam- 
bridge, UK, when CG and SJ visited in 2008, and at Institut Mittag-Leffler, Djursholm, 
Sweden, when SJ and AR visited in 2009. We thank Andreas Strombergsson for assis- 
tance with important references, and the referees for their helpful comments. 
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2 Notation, terminology and a summation theorem 



As mentioned above, G denotes a fixed connected multigraph with g vertices and no 
loops. For simplicity we assume that V(G) = [g] := {1, . . . , g}. We denote the number 
of edges in G by h. (Often we assume G to be <i-regular, and then h = dg/2.) Let 
A = Aq be the g x g adjacency matrix of G and let A = Aq be the incidence matrix of 
G, with g rows and h columns. Thus 

AA T = A + D G , (2.1) 

where Dq is the diagonal matrix with entries dc(i), i G V(G). Denote the eigenvalues 
of A by «i, . . . , ot g . 

In Section H] we also need a directed incidence matrix for G. Give each edge in G 
an (arbitrary) direction, and let Aq be the corresponding directed incidence matrix. In 
other words, Aq is the g x h matrix obtained from A by changing the sign of one of the 
two l's in each column. Then 

AqA t g =D g -A. (2.2) 

Our version of Laplace's method (Theorem 12.31) involves lattices. A lattice is a 
discrete subgroup of Mr. (Discrete means that the intersection with any bounded set 
in R N is finite.) It is well-known that every lattice £ is isomorphic (as a group) to 
Z r for some r with < r < n. The integer r is called the rank of £ and is denoted 
by rank(£). In other words, every lattice £ has a basis, i.e. a sequence Xi, . . . ,x r of 
elements of £ such that every element of £ has a unique representation Yll=i n i x i with 
rii G Z. Furthermore, the basis elements linearly independent (over R); 

thus the rank equals the dimension of the linear subspace spanned by £. 

The basis is not unique (except in the trivial case r = 0); if 5 = is any r x r 
integer matrix such that the determinant det(S) = ±1 (which is equivalent to the 
condition that both 5 and E" 1 are integer matrices) and ( is db basis of £, then 
Vi = J2j iij x j defines another basis yi, . . . , y r ; conversely, given (x«)i, every basis of £ is 
obtained in this way by some such matrix S. 

A unit cell of the lattice £ is the set {J2i=i^i x i : < < 1} for some basis 
of £. If £ C R^ has full rank N, and U is any unit cell of £, then {x + U} xe c is a 
partition of R^. 

The unit cells of a lattice £ all have the same r-dimensional volume (Hausdorff 
measure), where r = rank(£); this volume is the determinant (or covolume) of £, and 
is denoted by det(£). 

If (xi) r i=1 is a sequence of vectors in M. N , the symmetric matrix ((xj, Xj))l - =1 of their 
inner products is called their Gram matrix. It is well-known that linearly 
independent if and only if the Gram matrix is non-singular, i.e., if and only if the Gram 
determinant det((xi, x-,-))[ J=1 ^ 0. 

The following results are well-known. 

Lemma 2.1. If ( Xj^j^—-^ IS Qj basis of a lattice £ in M> N , then 

det((x 4 ,x J ))^ =1 = det(/:) 2 . (2.3) 
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Lemma 2.2. If Ci C £ 2 o^e two lattices of the same rank, then C2/C1 is a finite group 
of order det(£i)/ det(£ 2 )- 

The Hessian or second derivative D 2 <f>(xo) of a function at a point xq G R w is an 
N x N matrix; it is also naturally regarded as a bilinear form on R . In general, if B is 
a bilinear form on R^, it corresponds to the matrix {B{e^ . =1 , where (ej)^ is the 

standard basis. We define the determinant det(-B) as det (B(e i ,ej)) i j =v and note that 
if z\, . . . , zn is any basis in W N , then 



det (B(zi,Zj))f 
det({zi, Zj))^ =1 



(5) = "'\;""7;y =l . (2.4) 



We are interested in the restriction to a subspace. If B is a bilinear form on 1SL N 
and V C M. N is a subspace, we let det(B\v) denote the determinant of B regarded as a 
bilinear form on V. By (12.41) . this can be computed as 

det (B(zi, z i )) r . , 

for any basis zi, . . . , z r of V. 

We now state our general theorem for performing summation over a lattice using 
Laplace's method. 

Theorem 2.3. Suppose the following: 



(ii 
(iii 

(iv 

(v 

(vi 

(vii 

(viii 
fix 



C C M. N is a lattice with rank r < N. 

V C M. N is the r- dimensional subspace spanned by L. 

W = V + w is an affine subspace parallel to V, for some w G M. N . 

K C M. N is a compact convex set with non-empty interior K° . 

: K — > R is a continuous function and the restriction of <fi to K H W has a 
unique maximum at some point Xq G K° D W. 

<p is twice continuously differentiable in a neighbourhood of Xq and H := D 2 <P(xq) 
is its Hessian at x - 

ip : K\ — > R is a continuous function on some neighbourhood Ki C K of xq with 

Tp(x ) > 0. 

For each positive integer n there is a vector i n G R N with £ n /n G W, 

For each positive integer n there is a positive real number b n and a function 
a n : (£ + £ n ) H nK — > R snc/i i/iai, as — > oo, 

a n (£) = O (\e n<K£/n)+o(n) ) , £ G {C + 4) n ntf, (2.6) 
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and 



a n (£) = b n (^(£/n) + o{l))e n ^' n \ £ G (£ + £ n ) n n^, 

uniformly for £ in the indicated sets. 
Then provided det(—H\v) ^ 0, as n — ► oo, 

V flf!)~ (27T) r / 2 ^(x ) r/2 n*(*o) M 7 ) 

^ n[l) det(£)det(-H\ v y/* bnn 6 ' {ZJ} 

ee(c+e n )nnK v y v ,v ' 

We remark that Theorem 12.31 can be generalised to allow n to tend to infinity along 
any infinite subset I of the positive integers, with the same proof. (Then (viii) and (ix) 
need only hold for every n G I.) 



3 Expected number of perfect matchings 

A fractional perfect matching of the multigraph G is a function / : E(G) — > [0, 1] such 
that 

J2f( e ) = 1 for a11 v e v (°)- 

Note that every <i-regular multigraph has a trivial fractional perfect matching obtained 
by giving each edge weight 1/d. We often treat / as a vector (f(e)) ee E(G)- 

First, note that if there is a perfect matching at all in a lift L n (G) of G, then 
there exists a fractional perfect matching f of G such that nf(e) is an integer for each e. 
Indeed, suppose that M is a perfect matching of a lift of G. Let £ e be the number of edges 
from the fiber F e in M, for each edge e G E(G). Then the function / : E(G) — > [0, 1] 
defined by /(e) = £ e /n is a fractional perfect matching of G. Conversely, suppose that 
there exists a fractional perfect matching z = (z e ) e in G such that nz e is an integer for 
each e. We may construct an n-lift of G that contains a perfect matching as follows: 
First take nz e edges above each edge e G E(G), with all their endpoints disjoint. This 
yields n endpoints above each vertex i G G, so we have constructed the sets Vi, and 
a perfect matching. Extend this perfect matching to an n-lift by adding further edges 
between Vi and Vj for all edges e = ij. Consequently, L n (G) has a perfect matching 
with positive probability if and only if there exists a fractional perfect matching z with 
nz integer-valued. In the sequel, for a given graph G we consider only those values of n 
for which this holds, since otherwise trivially Xq = 0. 

Remark 3.1. It seems to be an interesting problem to characterize the set of such n 
for a given graph, but this is outside the scope of the present paper, and we note only 
the following examples: If G itself has a perfect matching then every n is allowed. On 
the other hand, if g is odd, then only even n are possible. If G is of odd order and 
hamiltonian, then the set of allowed n is exactly the set of positive even integers. If G is 
(i-regular, then (1/d, . . . , 1/d) is a fractional perfect matching, so every multiple of d is 
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an allowed n (but there might be others too). The result by Linial and Rozenman [12] 
implies that for a large class of graphs defined there, every large even n is allowed. Note 
finally that if n\ and 712 are allowed, then so is n\ + n-i- Hence the set of allowed n is 
always infinite, unless it is empty, so it makes sense to talk about asymptotic results. 

Suppose that there exists a fractional perfect matching z = (z e ) e in G with nz an 
integer vector. If a perfect matching in L n (G) has £ e edges in the fiber F e over e, then 

YleBv^e — n — n ^2 e ^v z e f° r every e, so (£ e ) e — nz belongs to the lattice C^q in M. E( - G ^ 
defined by 

£g) ;= |(i/ e ) e e Z E( - G) : = for every v G V(G)} 

= {ve Z E{G) : Au = 0}. 

(The superscript 1 denotes the first moment.) Here, and elsewhere when convenient, 
we think of the vectors as column vectors although we write them as row vectors for 
typographical reasons. Conversely, if £ = (£ e ) e is a vector such that £ — nz G Cq\ then 
£ is an integer vector and J2 eBv £ e = J2e3v nz e = n f° r ever y v - 

Given such an integer vector (£ e ) e G Cq + nz, let us compute the expected number 
of perfect matchings in L n (G) with £ e edges in the fiber F e . Clearly this number is zero 
unless < £ e < n for all e. Then the endpoints of the edges in the matching may be 
chosen in 

ways, and for each choice, there are £ e \(n — £ e )\ possibilities for the fiber F e , with 
probability \jn\ each. Hence, defining K = [0, 1] E ^ we have 

E(X G )= Yl a ™W t 3 - 1 ) 

ie(C^ +nz)r\nK 

where 

(Recall that h denotes the number of edges in G.) 

We wish to evaluate the sum (13. ip asymptotically by Laplace's method: more pre- 
cisely, by applying Theorem 12.31 We use Stirling's formula in the following form, valid 
for all n > 0, where x V y :— max(a;, y), 

ln(n!) =nlnn-n+iln(nVl) + iln27r + 0(l/(n + l)). (3.2) 
Let x e = £ e /n for all e G E(G). Applying (13. 2p we obtain, uniformly for £ G (Cq + 
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nz) fl nK, 

Ha n (£)) = (g-h)\n(n\)+ ]T (ln((n - 4)0 - WJ) 

e€E(G) 

= (g-h) (n(ln(n) - 1) + \ ln(n) + | 1h(2tt) + 0(1/ n)) + E ( n ~ 2 4)(hi(n) - 1) 

e€E(G) 

+ n ^2 ((1 - £ e ) ln(l - x e ) - x e ln(x e )J 

e£E(G) 

+5 E N(i- I ,)v n -')-^v n -' ) )+ E (4TT + ^dm) • 



eeS(G) ee-B(G) 

SlIiCG 

E '. = §E£'« = IE» = k»> 

eGE(G) V e^v v 

after cancellation, a n (£) can be expressed as 

1 ) un £ e + 1 / ^ ( n — max 4 I < 
where, for a; G R E< - G - ) , 

6 n := (2t™)^-^/ 2 , (3.3) 
0(x) := E(( 1 - x e) ln ( 1 - x e) -x e \n(x e )), (3.4) 

e 

1 — 2; X 



a n (£) = b n iP(£/n) exp(ri(f)(£/n)) (l + o( . ) ) + O 



^)-n(^fJ ' (3 - 5) 

except that if some x e or 1 — x e is 0, we replace it by 1/n in (13.51) . This implies that 
a n (£) satisfies condition (I2.6P of Theorem 12.31 with the above b n , 4>, and ip. We will now 
check all the remaining assumptions of Theorem 12.31 Let 

W := jx = (x e ) G M E(G) : J^x e = 1 for every u G = {x : Ax = (1, . . . , 1)}. 



e3v 



As is well-known, and described in Section [H] in detail, the sum (13.11) is dominated by the 
terms where <p(£/n) is close to its maximum. In order to find the maximum, we restrict 
ourselves to regular multigraphs, where the result is simple. (The method applies to 
other graphs as well, provided one can find the maximum point (s) of 0.) 

Lemma 3.2. Suppose that G is d-regular, where d > 3. Then <p defined by (13.41) 
has a unique maximum on K fl iy={x G K : Ax = (1, . . . , 1)}, attained at the point 
x° = (1/d, . . . ,1/d). The maximum value is 

m, ox 9, {(d-iy- 1 
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and, for if) in (13.51) and the Hessian D 2 cf), 



^) = {d-lf'\ D^(x°) = -^-^-I. 
Proof. We write <p = \ J2 ve v(G) <^> where 

<f) v (x e :e3u) = ^((l - x e ) ln(l - x e ) - x e ln(x e )). (3.6) 

Fix a vertex v G V(G). We rename the variables x e , e 3 v, by xi, . . . , Xd, for convenience. 
Since <ft v is continuous, it has a maximum over the compact set 



a 

£ d :={(£,), G[0,l] d :5> l = l}. 



Let x v 6 S,j be a maximum point of <p v . Assume first that x v is an interior point, i.e., 
that x v G (0, l) d . Then the function f(y) = 4> v (x\ + y,x v 2 — y,x%, . . . ,Xj) achieves a 
maximum at y = 0. Therefore, /'(0) = and by the chain rule, 

dMx± _ g^c) 

By the same argument (or by the general Lagrange multiplier method), we have that 
for some constant c v > 

-^^—( xV ) = c v , for i = 1, . . . , d. 

But 

— (x ) = - ln(l - Xi) -Inxi- 2, 

OXi 

so 

2^(1 — xi) = exp{— c v — 2} for alH = 1, . . . , d. 

This implies that the x^s are all at the same distance from 1/2. That is, for some 
constant c' v > we have x\ = 1/2 ± c' v for i = 1, . . . , d. Since ^ x\ — \ and ti > 3, 
we have to choose the minus sign for all i, and thus all x\ are equal. Since x v G E^ we 
conclude that x\ — 1/d for i — 1, . . . , d. 

We also have to consider the boundary of E d . If, say, x\ = and < x| < 1, 
then / above is defined for small positive y with /'(0+) = +oo, so x v cannot be a 
maximum point on E^. The only remaining points are those with all Xi G {0, 1}, but 
then 4> v (x) = 0, while <f> v (l/d, . . . , 1/d) > 0, so these too cannot be (global) maximum 
points. Hence x v is the unique maximum point for (p v on E^. 

Setting x° = (1/d,..., 1/d) G W, we have for all x G K n W, 

4>{x)<\Y,Ux v )=<t>{A- 
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Moreover, the inequality is strict for all x ^ x°. This proves that x° is a unique maximum 
point of in K D W . Clearly, x° belongs to the interior of K. Moreover, and 
i/j(x°) are given by the formulas stated in Lemma [3.21 

Finally, the Hessian D 2 <ft(x) is diagonal with entries (1 — x e ) _1 — x~ . Hence, at x° 
we have D 2 (j)(x°) = -^f^- □ 

We have verified all assumptions of Theorem 12.31 f° r an Y neigbourhood K\ of x° 
with K\ C K° . To apply formula (12.71) . we still need to compute the rank of the lattice 
Cq and its determinant det(£^). 

Lemma 3.3. (i) // G is non-bipartite then the lattice C^q has rank h — g and deter- 
minant det(£g } ) = | det(A + D G ) l l 2 . 

(ii) If G is bipartite then the lattice C^q has rank h—g+l and determinant det(£^' ) ) = 
det(A' + D' G ) l l 2 , where the matrix A' (respectively, D' G ) is obtained by deleting the 
last row and column of A (respectively, Dq). 

Proof. For v G V{G) define the vector x v = (l[v G e], e G E(G)\ given by the row 
of the incidence matrix A corresponding to v. For convenience, rename these vectors 
X\, . . . ,x g . Then, by (12. ip . the Gram matrix of xi, . . . ,x g is AA T = A + Dq. This 
matrix is singular if and only if there exists a non-zero vector y — (y v ) G R^( G ) with 
yA = 0. This is equivalent to yi = —yj for every edge ij, and it is easily seen that, when 
G is connected, such a non-zero vector y exists only if G is bipartite, and that if G is 
connected and bipartite, there is a one- dimensional space of such solutions y. 

Consequently, in the non-bipartite case (i), the vectors linearly inde- 

pendent. We apply Lemma 16.21 with N = h, m = g and using the vectors x%, . . . ,x g . 
Let C, C L and Co be as in Lemma [6.21 Then C^q = C 1 , and thus Cq has rank h — g, 
by Lemma [6.21 Furthermore, by Lemma [2. II and (12. ip . 

detOC) = (det((x J ,x,))f J=1 ) 1/2 = det(A + D G ) 1 ' 2 . 

Moreover, (t v , v G V(G)) solves (16.11) if and only if t v = —t w (mod 1) for every edge vw. 
Going around an odd cycle, we see that t v = or t v = 1/2 for every vertex on the cycle. 
Since G is connected, it follows that there are exactly two solutions to (16. ip : t v = for 
every v and t v = 1/2 for every v. Hence q = 2 in Lemma [6 .2\ and the result follows. 

Now suppose that G is bipartite. Then the vectors linearly inde- 

pendent and be written as a {±l}-combination of X\, . . . since the sum 

of vectors x v over all vertices v on either side of the vertex bipartition gives the vec- 
tor (1, 1, . . . , 1). We apply Lemma 16.21 with N = h, m = g — 1 and using the vectors 
xx, . . . , x g -\. The lemma asserts that Cq = C 1 - has rank h — g + l, and 

det(£ ) = (det((x l ,x J ))f-l 1 ) 1/2 = det(A' + D ! G fl 2 . 

Finally, let w G V{G) correspond to x g . If (t v , v G V(G) \ {w}) solves (16. ip then t u = 
for every neighbour u of w. In turn this implies that t u = for every vertex u at distance 
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2 from w, and iterating this shows that t u = for all vertices u in the connected graph 
G. Therefore q = 1 in Lemma [6.21 and the proof is complete. □ 



Example 3.4. When G = K 4 , 

det(A + D G ) 

Thus Lemma I3.3( i) says that C^q has rank 2 and 

(ih 



3 111 
13 11 
113 1 
1113 



48. 



detect) 



12. 



Example 3.5. Let G = K\ be the multigraph with two vertices and three parallel 
edges. Then A + Dq = I J and deleting one row and column gives the lxl matrix 

(3). Hence Cq has rank 2 and det(£^) = \/3, using Lemma ETBT ii). 

We are ready to apply formula (12.71) of Theorem 12.31 

Theorem 3.6. Suppose that G is d-regular, where d > 3. 

(i) If G is non-bipartite then 

2 (d- 1)*/* ( d-\ \ dg/4 ~ 9/2 /(d-l) d - 1 \ W2 



EX 



G 



Vdet^ + dT) Vd(d-2) 

2 (d_ i)(d-i)g/2 ( {d-\)' 
(d(d - 2)) rf 3/4- 9 /2 v /det(A + ^fy V d d ~ 2 



d d " 2 



(ii) // G zs bipartite then 

(d - l) d °/ 4 



EX 



G 



d- 1 



y/detfA' + dl) \d(d-2 

(d- l)(d-l)9/2+l/2 



dg/4-g/2+l/2 



: (2vrn) 



(27m) 1 / 2 
1/2 ( (d 



(d 



cZ-l\ W2 



1 



d rf " 2 
\<2-l\ W 2 



i-2 



(d(d - 2))^/4- 9 /2+i/2 ^det^' + d/) 
where A' is obtained by deleting the last row and column of A. 

Proof. Let r be the rank of £q\ and recall that the Hessian H = D 2 (f)(x°) is diagonal and 
equals — d ^Zi I by Lemma [3T2l Thus H\ v = — d ^~^ I too, and det(— H\v) — { d ^d~? Y ■ 
Hence the result follows from fl3.ll) and Theorem 12.31 using Lemmas 13.21 and 13.31 and 
the fact that h = dg/2. □ 
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Example 3.7. For G = K±, d = 3, g = 4 and thus, using Example 13.41 



4 / i \ 2n / j \ 2n 



EA' G 2 ' 2 /4V 8 /4 



3^ V 3 / 3\/3 V 3 



Example 3.8. For the bipartite multigraph G = K% with two vertices and three parallel 
edges we have d = 3, g = 2 and by Example I3.5[ 

8 , /4 



E Xr; ~ — 7= v vrn - 
3v/3 V 3 

4 The second moment of Xq 

We now work towards an asymptotic expression for the second moment of Xq, using 
the same approach as in the previous section. To simplify our calculations we consider 
only regular multigraphs G of degree at least three. 

Given a pair (Mi, M 2 ) of perfect matchings in L n (G), for a vertex i £ V(G) and two 
(possibly equal) edges e, / 3 i, let £j e / be the number of vertices in Vi whose incident 
edges in Mi and M 2 lie, respectively, in the fibers F e and Ff. Form these numbers into 
the gd 2 - dimensional vector £ = £(M\, M 2 ) = (ii e f : i £ [g], e, f 3 i). Let 

V* := | (zi e f : i £ [g], e, f 3 i) £ M 9a!2 : for every e £ E(G) with endpoints z and j, 

Ziee — z jee, 2 «e/ = ^je/; Z i/e = Z .j/eJ- 

Then the vector £ belongs to the set 

Q := {(z ief ) eV*H Z 9d2 : J2 z ^f = n for IE [<?]}. 

(The three conditions in V* follow from consideration of the edges in Mi D M 2 , Mi and 
M 2 , respectively.) Fix a particular vector z with nz £ Q. (By our assumption that there 
is a perfect matching in L n (G), it follows that at least one such vector exists.) Then 
Q = + nz, where C^) is the lattice defined by 

4 2) := {(vief) U> S : v *f = for z £ [<?]}. 

(The superscript 2 denotes the second moment.) 

Given a pair (Mi,M 2 ) of perfect matchings and thus a vector I £ Q, we further 
define, for an edge e £ -E(G) and an endpoint i of e, 

>5e ^ie(^) ^ ^ ^iefj t e ^je(^) ^ ^ ^ifei ^ie(^) ^ ^ ^iff'i 
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these are the numbers of edges in the fiber F e that belong to M\ \ M2, M 2 \ Mi and 
(Mi U M 2 ) c , respectively, so they do not depend on the choice of endpoint i of e. We 
have, for every edge e and endpoint i, 

S e ~\~ t e ~\~ U e ~\~ £iee 

We now calculate the expected number of pairs of perfect matchings (M l5 M 2 ) in 

(2) 

L n (G) corresponding to a given nonnegative integer vector t = (£j e /) G C G +nz. First, 
partition each Vi into d 2 subsets of sizes (A'e/)e,/9ij this can be done in 

i=1 lle,/3iW- i=le,/9i 

ways. Given these partitions there are 

S e \ t e \ U e \ £iee- 

possibilities for the fiber F e (where i is an endpoint of e), with probability 1/n! each. 

Hence the expected number of pairs (M\,M2) of perfect matchings in L n (G) which 
correspond to the vector I is given by 



„w= n ^n(n( f T^) 1/2 n A 



Thus we can write 

££(£^ 2) +n 2 )nn_ft' 

where K = [0, 1] 96 ' 2 . This will allow us to apply the same arguments as used in Section [31 
We now switch to continuous variables x G M sd2 , where corresponds to £i e f/n. 
Define the functions <Ji e = <Ji e (x), r ie = r ie (x) and 7j e = ji e (x) to be continuous scaled 
analogues of Sj e , t ie and Ui e respectively. That is, 

(Tie ^ %iefi Tie ^ %ifei ''lie ^ •^ijf'i 

flhf^e fBiJ^e f,f'Bi;f,f'^e 

so that (iie(£/n) = Si e (£)/n and so on. Then, applying (13.21) . it follows that a n (£) satisfies 
condition (12.61) of Theorem 12.31 with 



f>„ = (2nn)" 2+3, " 2 - d '«l 2 , 



*w=nn(^) 1/4 11 <r. 



X ' 

i£[g] eBi x lco / /9i,/^ e 



(f)(x) 



ie[g] e3 



2 XieflnXiefY (4.2) 
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(Again, if some Xj e /, <T ie , Tj e or 7j e is 0, then we replace it by 1/n in the definition of 
Let be the domain defined by 

W := {(Xief) eV* x ^f = 1 for i e Id}}- 

e,f3i 

We conjecture that for all connected d-regular multigraphs G with no loops, the function 
has a unique maximum on K H VF, attained at the point 

x° = (l/d 2 ,...,l/d 2 ). 

Unfortunately, we have been unable to prove this, and have only been able to verify this 
computationally for d = 3. For future reference, note that 

= ((d - l)^ 2 f , <K*°) = In ( ^r^ ) • (4-3) 

One approach to finding the maximum of is to mimic the proof of Lemma I3~^l The 
function <fi can be written as the sum over i — 1, . . . , g of functions 0j, where the sets of 
variables appearing in different 0; are disjoint. For convenience we drop the index i and 
rename all variables corresponding to vertex i as x e f := Xi e f, and let a e := a ie , r e := r ie , 
7 e := 7 ie . Then 

= | -j cr e In cr e + r e In r e + 7 e In 7 e -x ee In x ee - 2 ^ x ef \nx ef j-. 

Since G is d-regular and 0j depends only on the degree of i in G, all the functions 0j 
are equivalent under relabelling of variables. 
Now define the domain 



X<P = {K/)e,/3i e [0, if : 53 x e/ = l| 



It suffices to prove that 0j has a unique maximum on E,p attained at the point (1/d 2 , . . . , 1/ d 2 ) 
Applying the Lagrange multiplier method to S rf 2, we see that at an interior maximum 
point, all partial derivatives of 0j must be equal. This gives d 2 — 1 (non-linear) equations 
(together with ^ e . x e f = 1) to be solved for d 2 variables. We tried to solve this system 
using Maple. Unfortunately, Maple seems unable to handle the computations for d > 4. 
Hence we only have the desired result for d = 3. 



Lemma 4.1. If G is 3-regular then the function <fi defined by ( \4-2 ) has a unique maxi- 
mum on K C\W attained at the point (1/9, . . . , 1/9) G M. 9g . 

Proof. As explained above, we consider only the function 0, for a fixed vertex i. Using 
Maple, we solved for points in { (x e f) e j : J2 e / x ef = l} where all the 9 partial derivatives 
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of <f>i are equal. Exactly four solutions were found, of which only one lies in [0, l] 9 , giving 
the point x° = (1/9, . . . , 1/9) G E 9 . (The other three solutions each contain both positive 
and negative entries.) We have <j)(x°) = ln(4/3). 

It remains to consider the boundary, where one or several x e f = 0. If x ee = and 
7/ > for / 7^ e, then -^—(j)(x) = +oo, and thus x is not a maximum point. Similarly, 
x cannot be a maximum point if x e f = 0, where / and at most one of a e , Tf and 7/' 
(where /' is the third index) vanishes. It is easily seen that the only remaining cases are 
when the only non-zero variables (after relabelling the indices as 1,2,3 in some order) 
are {212,^21}, {^11,^22,^33} or {^11,^12,^13}, or a subset of one of these. In the first 
case we have = 0. In the two latter cases, <pi equals, after relabelling, |0„ defined in 
(13.61) (at the corresponding step of the first moment calculation), and thus the maximum 
over one of these sets is |ln(4/3) < 0(xo). (We omit the details.) Hence, there is no 
global maximum on the boundary. 

Consequently, x° is the unique maximum point of 0j on Eg. Arguing as in Lemma I3T21 
completes the proof. □ 

Let V = W — z be the subspace spanned by C G , i.e., 

V := {(z fe /) e V* : x ^f = for i E [#]}. 

Theorem 4.2. Suppose that G is d-regular, where d > 3. If the function <fi defined in 
(14.21) has a unique maximum on K C\W at x° = (1/d 2 , . . . , 1/d 2 ), then 

E (X 2 G ) ( ( f ~ l)dd ~ 2 ^ (27m )^+^+^/4- rf V2 ( V-Jf 1 )'", 

V G) det(4 2) )det(-tf| y )V2 l ^ V J 

where r is the rank of C^q and H = D 2 (f)(x°) is the Hessian of <f> at x°, provided the 
determinant in the denominator is non-zero. In particular, this expression holds for all 
3-regular connected graphs G. 

Proof. This is now an immediate consequence of Theorem I2.3[ using (14. ip and (14. 3 p . 
The final statement follows from Lemma [4.11 □ 

(2) 

It remains to calculate the determinants of C G and —H\v, and the rank r. In the 
non-bipartite case, part of this is covered by the next lemma. 

Lemma 4.3. Suppose that G is non-bipartite and d-regular, where d > 3. Recall that 
h denotes the number of edges in G, so h = dg/2. Then the lattice C G has rank 
d 2 g — (g + 3/i) = d 2 g — g — 3dg/2 and determinant 

det(4 2) ) = 2* h ' 2 -^/ 2 - 2 (d{d - 2)) hl2 - g/2 det(dJ + A) det(d(2d - 3)1 - A) 1 ' 2 

= 2 3h / 2 - 3 ^ 2 ' 2 (d(d - 2)) h/2 - g/2 f[(d + ai )(d(2d - 3) - a,) 1 / 2 , 

i=l 

where a±, . . . ,a g are the eigenvalues of A. 
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Proof. The linear space V spanned by c!q is the subspace of M. 9d2 orthogonal to the 
following g + 3h vectors: 

• one vector x°i for every j G V(G), with xf e t = l[i = j\. 

• one vector x le for every e G E(G), with x\lj = a i£ l[e = f = e}. 

• one vector x 2e for every e G E(G), with = aj e l[e = e^/]. 

• one vector x 3e for every e G E(G), with sfjy = a i£ l[e ^ e = /]. 

Relabel these vectors (in this order) as Xi, . . . , x 9 +3/i- Then their Gram matrix T can be 
written in block form, with blocks of dimensions g,h,h, h: 



( d 2 I 
A T 
(d-l)A T 
\(d~l)A T 



A (d-l)A 

21 

2(d-l)7 

A T A-2I 



(d-l)l\ 





A T A-2I 
2(d- 1)1 



In order to evaluate the Gram determinant det(T), we may make an orthogonal change 
of basis in the first component M 9 , and another orthogonal change of basis in each of 
the components M. h (we choose the same change in all three). It is well-known that we 
can make such changes of basis such that any given g x h matrix B obtains the form of 
a diagonal g x g matrix D s with h — g additional columns of O's; this is known as the 
singular value decomposition of B, and is easily seen by choosing an orthonormal basis 
Zi, . . . , Zh itl M. h such that B T B is diagonal, and then choosing an orthonormal basis in 
M 9 containing the vectors Bzi/\\Bzi\\, for all i such that Bzi ^ 0. We choose such bases 
for B = A. The diagonal entries si . . . , s g of D s can be assumed to be non-negative, 
and they are identified by the fact that the eigenvalues of BB T = AA T are {s 2 }. By 
( 12.21) . we thus have 

s 2 = d- a;. (4.4) 
(D s , 0) a g x h matrix with non-zero elements given by (14.41) . 



Hence, with D s 



detr 



d 2 I D s (d-l)D s (d-l)D s 

Dj 21 

{d-l)D T s 2{d-l)I D T S D S -2I 

(d-l)Dj D T s D s -2l 2{d-\)I 



(4.5) 



Since D s is a diagonal matrix, we can reorder the rows and columns in (14. 5 p so that we 
obtain a block diagonal matrix with g 4 x 4 blocks 



/ d 2 Si (d-l)si (d-l)si\ 

Si 2 

{d-l)si 2(d-l) sf-2 

\{d-l)si s 2 -2 2{d-l)) 



(4.6) 
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and h — g identical 3x3 blocks 



/2 \ 

r := 2(d-l) -2 . (4.7) 
\0 -2 2(d-l)/ 

Hence, by straightforward calculations, 

9 

det(r) = det(r )^ 9 JJdet(r 4 ) 
1=1 

= (8rf(rf-2))^] £ [(2rf- Sl 2 ) 2 (2d 2 -4rf + S , 2 ) (4.8) 
s 

= (8d(d - 2)) h ~ 9 \\(d + a l ) 2 {d(2d - 3) - on) 

8=1 

Since G is non-bipartite, — d < a>i < d for every i, and thus ( 14. 8 p shows that det(r) ^ 0. 
Hence, the vectors x\, . . . , £ g +3h, or in different notation 

{x 0j : j E V(G)} U {x l£ , x 2e , x 3£ :ee E(G)}, (4.9) 

are linearly independent, so they form a basis in V . 

We apply Lemma \6.2\ with iV = d 2 g, m = g + 3h = g + 3dg/2, and using the vectors 
xi,..., x g+3h in flMD. Then C {2) = C L . Hence, rank(£[j 2) ) = N — m = d 2 g — g — 3h. We 
have det(£o) = det(r) 1 ' 2 by Lemma 12.11 Finally, we claim that there are 4 solutions 
(mod 1) to (16.11) : if we let t j denote the coefficient of x - 7 , and so on, the solutions have 
toj = to for all j and t\ e = ti, t 2e = t 2 , t 3e = t 3 for all e, where (to, ti, t 2 , fa) = (0, 0, 0, 0), 
(0, 0, |, |), (|, |, -, 0), or (|, |, 0, |). (To prove this, first consider the equations in ( 16.11) 
which correspond to variables Xi ee , and use the existence of an odd cycle. This gives the 
possible values of to and t\. The rest of the proof follows by considering the equations 
in (16.11) which correspond to variables Xi e f for a given vertex i, with e ^ /.) 

Hence q — 4, and Lemma 16.21 yields 

det(£g ) ) = det(£ ± ) = det(r) 1/2 /4. 
The result follows by lOj) . □ 

Example 4.4. For G = K4, we have d = 3, g = 4, h = 6, and A has the eigenvalues 
3,-1,-1,-1. Hence Lemma yields det^) = 2 7 3 5 / 2 5 3 / 2 . 

We believe that there is a similar result for regular bipartite graphs, but we have not 
explored it. (Presumably, the rank is then d 2 g — g — 3h + 2). 

Unfortunately, we have not been able to find a similar general formula for det(— H\y) 
in Theorem 14.21 However, this quantity can be calculated directly for a particular graph 
G, once a basis for Cq is known. 
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Example 4.5. When G = K±, using Maple we found a basis {z±, . . . , Zu} of V and 
then calculated det(-H\ v ) = 2~ 22 3 28 5 _1 ll 3 using (|2~5|) . Hence by Theorem [O and 
Example 14.41 

E(X 2 )~ 2 16 3- 9 / 2 5- 1 ir 3 / 2 

Example 4.6. When G = K\ is the multigraph with two vertices and three parallel 
edges, Maple computations confirmed that L^) has rank 9 and gave det(£^) = 2 4 3 3 / 2 
and det(-H\ v ) = 2~ 16 3 18 5 2 . Hence by Theorem S^l 

E(X 2 )~2 11 3- 9 / 2 5- 1 t™ 

5 Short cycles in random lifts 

Let Zk denote the number of cycles of length k in L n (G), for k > 2. (Note that Z 2 is zero 
unless there are multiple edges in G.) To apply the small subgraph conditioning method 
to Xq-i we must understand the distribution of short cycles in random lifts, as well as 
their interaction with perfect matchings. This will enable us to verify conditions (Al) - 
(A3) of [HI Theorem 9.12], with their Y n given by our X G (the index n is suppressed), 
and with their X^ n given by our Z^. 

To compute the limiting distributions in (Al) and (A2) of [11, Theorem 9.12], we 
will use the method of moments. Moreover, for (A2) we will be guided by [TTJ Lemma 
9.17 and Remark 9.18], which tell us that we need only compute asymptotically 

E(X G (Z 2 ) j2 ---(Z m ) jm )/EX G , 

for integer constants m > and j 2 , ■ ■ ■ ,j m > 0. Here (Z)j denotes the falling factorial 
Z(Z-l)-..(Z-j + l). 

Let k be a fixed positive integer. It is more convenient to count rooted oriented 
fc-cycles, which introduces a factor of 2k into the calculations. A /c-cycle in L n (G) 
can be then thought of as a lift of a non-backtracking closed A;- walk in G, which is a 
walk ioeiiie 2 ■ ■ ■ %k-\^k in G such that tj is an edge of G with endpoints {ij,ij + i} and 
Cj 7^ ej-i, for 1 < j ' < k. (Here and throughout this section, arithmetic on indices in 
k- walks is performed modulo k.) Note that if G is simple then any three consecutive 
vertices on the walk must all be distinct. These walks arise in various contexts (see for 
example [HO [10]) and have also been called irreducible [9] and non-backs cattering [13] . 
Denote by Wk the number of non-backtracking closed k- walks in G, for k > 2. 

The following lemma shows that condition (Al) of [HI Theorem 9.12] holds. 

Lemma 5.1. Let = Wk/(2k) for all k > 2, where Wk is the number of non- 
backtracking closed k-walks in G. Then Z^ ~ Po(Afc), jointly for all k > 2. 

Proof. Fix a non-backtracking closed fc-walk C = ioe\i\ ■ ■ -ik-iCk in G. The (oriented) 
A;-cycle C' = f\f 2 ■ ■ ■ fk in L n (G) is a lift of C if fj G F ej for j = 1, . . . , k. Hence the 
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number of possible lifts C of C is (1 + o(l))n k , and each will appear in L n (G) with 
probability (1 + o{l))n~ k . It follows that 

E z * = E E p ( c " c = 1 + 

C C" 

Similar arguments hold for higher joint factorial moments, completing the proof. □ 

For the remainder of this section we restrict our attention to d-regular multigraphs 
with d > 3. Next we verify condition (A2) of pTTl Theorem 9.12] using the approach 
suggested in [HI Remark 9.18]. 

Lemma 5.2. Suppose that G is d-regular with d > 3, and for k > 2, let 

w= ( 1 + (i^r) J A " 

Then for any integer m > 2 and non-negative integers • • • ,jm, 
E(X G (Z 2 ) j2 ■ ■ ■ (Z m ) jm ) 



EX 



G 



i=2 



J% as n —>■ oo. 



Proof. For ease of notation, throughout this proof we write P(M) := P(M C L n (G)), 
P(M,C") := P(M C L n (G),C C L n (G)), and so on. First we estimate E(X G Z fc ). We 
write 

e(x g ^) = EEE p ( M ' c ') = E p ( M ) E E nc\M), 

m c a m c a 

where the sums extend over all possible perfect matchings M in L n (G), all non-backtracking 
closed fc-walks C in G, and all their possible lifts C, respectively. 

To calculate the inner double sum, we fix a perfect matching Mo and condition on its 
presence in L n {G). Let C = ioeiii . . .ik-iek be a given non-backtracking closed fc-walk 
in G. For a lift C of C with edges /1/2 • • ■ fk, let 

(c0= fiifA 6 «b, tal<j<fc 

10 otherwise, 

To estimate the expected number of lifts of C given Mo, we break the sum over all 
C according to the vector £(C): 

j2nc\M )= E E nc\M ). 

C «e{0,l} fc C':£(C')=u 

Let £ e be the number of edges of M in the fiber F e , and say that M is good if 

— n/d\ < n 2 / 3 for every e. 
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We may assume that Mq is good, since the calculations for the expectation in Section [3] 
show that the contribution from other matchings is negligible. (Specifically, this follows 
from the proof of Lemma [6731 in particular the fact that S2 = o(l), S3 = o(l), using 
notation from that proof.) 

Hence, for a given u = (ui, u 2 , . . . , Uk) G {0, l} fc , 

¥(C'\M ) 



n — n/d 

Let tooiu) and toi(w) be the numbers of substrings 00 and 01 in u, respectively. Next we 
prove that the number of lifts C = f± • ■ ■ fk of C such that £ (C) — u is asymptotically 
equal to 

2n\ too(u} /n\*oi(«) 



Indeed, let V ie be the set of endpoints in Vi of the £ e edges in M D F e , for i incident to 
e G E(G). If, say, U\ = u 2 = 0, which means that both, /1 and f'2, are not in M , then we 
can choose the end of /1 in V ix from Vi 1 \(Vi iei UVi iea ), and |Vf 1 \(T^ iei UVi ie2 )| ~ n — 2n/d 
since we assume that Mq is good. Similarly, if m = and M2 = 1, which means that 
/1 ^ M but /2 G Mo, then we have to choose the end of /1 from Vi ie2 , a set of size 
~ n/d. Note also that if ui = 1 then we must have M2 = 0, and if we have already 
selected the end w of fi in Vi , then the other end of f\ is completely determined as the 
partner of w in Mq. 

Multiplying these two expressions together yields that 

J2 P(C"|M ) = b UlU2 ■ ■ ■ b Uk _ lUk b UkUl + o(l), 
where b 00 , 6 i, b w , fen form the matrix 
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d-2 1 
d-1 d-1 







Note that -B has eigenvalues 1 and —l/(d— 1). Summing over all u = (u\, . . . , Uk), we 
find that the conditional expected number of lifts of C is 



E = Tr(5 fc ) + o(l) = 1 + (— l -) 



k 

+ 0(1). 



Hence the expected number of fc-cycles in L n (G), conditioned on the existence of a given 
good perfect matching M , is asymptotically equal to 



«•/.• / , / -1 N ' 



2A; \ Vd-1 
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Finally, 

M 

All the above calculations work similarly for higher factorial moments and yield the 
desired result. □ 

Denote a directed edge of G by (e,i,j), where e G E(G) is incident to i,j G V(G) 
and i 7^ j; this denotes e directed from i to j. Now let i? be the dg x c?<? matrix with 
rows and columns indexed by directed edges of G, and 



^(e,ij'),(/,p,g) 



1 if p = j and / 7^ e, 
otherwise. 



(Here i? is the adjacency matrix of a version of the directed line graph of G, where 
[/-turns are forbidden.) Then 

w k = Tr(R k ) = 6 k + --- + 6 k dg} (5.1) 

where 81, ... , Q dg are the eigenvalues of R. Note that d — 1 is an eigenvalue of R with 
eigenvector (1,1,..., 1) T ; since R has non-negative entries, this is the eigenvalue with 
largest modulus. Now for k > 2, the quantity fi k defined in Lemma [531 equals 

■1 ^ k 



fi k = (1 + 5 k )Xk, where 5 k = [j—[J > -1- 
Therefore the quantity Afc5j? in condition (A3) of [TT], Theorem 9.12] is 

E Afc<** = fc (d-l) 2k = ^ ^ £ Ud-l) 2 ) 

k fc>l v 7 fc>l t=l vv ' ' 



2£f V 



which is finite as required. Furthermore, 
exp 



(E = ( d - l ) da (ri« d - x ) 2 - ^ 

= (d-l) d9 det((d-l) 2 I - R)~ 1/2 . (5.2) 

In order to assist with the verification of condition (A4) from from [11, Theorem 9.12], 
we will rewrite this expression in terms of the adjacency matrix A of G. The following 
result was proved by Friedman [9]. 
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Lemma 5.3. [9, Theorem 10.3] Suppose that G is d-regular with d > 3 and let 
ax, . . . ,a g be the eigenvalues of the adjacency matrix of G. For i = 1, . . . , g denote the 
roots of the quadratic x 2 — c*iX + d — 1 = by {if and j3~ . That is, 



A + = hi + v K - ( d - !)> A" = \<*i - v K - ( rf - !)■ 

T/ien i/ie eigenvalues of R are (3f, (3~ for i = 1, . . . , g, together with 1 and —1, the latter 
two repeated g(d — 2)/2 times each. Hence, for k > 2, the number of non-backtracking 
closed k-walks in G is given by 



\g{d - 2) (1 + (-l) fc ) + ((A + ) fe + Wt) ■ 



i=i 



Note that there may be repetitions among and some of these may coincide 

with ±1. Hence the multiplicities of these eigenvalues may not be exactly 1 or g(d—2)/2: 
see Example 15.51 below. 

We now use Lemma I5~3l to rewrite (15.21) in terms of the eigenvalues of the adjacency 
matrix of G. 

Corollary 5.4. Suppose that G is d-regular, with d > 3. The expression in (15. ,21) can 

be written as 

exp X k 5 l 

= (d - l) d9 -^ 2 ((d - l) 4 - l)-( d - 2 )s/ 4 det((d - l) 3 + 1)/ - {d - l)A)~ 1 ' 2 

9 

= {d- l) d9 - 9 ' 2 {{d - l) 4 - l)-( d ~%/ 4 Y[ {(d - l) 3 + 1 - (d - i)aiy 1/2 . 

1=1 

Proof. It follows from Lemma 15.31 that the characteristic polynomial of R is given by 

dg 9 
det(AI -R) = - 6i) = (A - l)( d - 2 ^/ 2 (A + lf d - 2 ) 9 / 2 J](A - /3+)(A - (3~) 

i=l i=l 

= {X 2_ 1) (d-2)9/2f[(\ 2 -a i \ + d-l) 
i=l 

= (A 2 - l)( rf - 2 )s/ 2 det((A 2 + d-l)I- XA). 
The proof is completed by substituting this into (15.21) with A = (d — l) 2 . □ 

Example 5.5. When G = K4 the eigenvalues of A are a± = 3, «2 = «3 = 0:4 = — 1. By 
Lemma |5~3| the eigenvalues of R are 2, 1 (three times), —1 (twice), and |(— 1 ± y/7i) 
(three times each), so the number of non-backtracking closed A;- walks in is 

^ + 3 + 2(-l)' + 3(ZI±^) 1 + 3 (ZI_^)' 
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Furthermore, by Corollary 15.41 



exp 



Afe^fe = 2 10 15- 1 det(9J - 2A)- 1 ' 2 = 2 10 3^ 5" 1 ll^ 2 . 



Example 5.6. The multigraph with two vertices connected by d parallel edges has 
adjacency matrix 

We have P^,P^ = ±(d — 1), ±1 and by Lemma 15.31 the matrix R has eigenvalues 
±(d — 1) and ±1, the latter with mulitiplicities d — 1. Hence Wk = 2(d — l) k + 2(d — 1) 
if k > 2 is even, and Wk = if k is odd. Corollary 15.41 yields, after some algebra, 



exp X k 6 l^J 



(d - lf^d-^id - 2)- d ' 2 {d 2 -2d + 2)- d ' 2+1 ' 2 . 



For example, when d = 3 this is 2 5 3 3/2 5 1 , while for d = 4 it is 2 15 / 2 3 7 5 3/2 . 

To complete this section, we prove a concentration result for the number of perfect 
matchings in L n (G) when G = K4 and when G is the multigraph K\ with 2 vertices 
and 3 parallel edges. We conjecture that the analogous result is true for any connected 
(i-regular multigraph G with no loops, where d > 3, with 5k = — (l/(d — l)) fc . 

Corollary 5.7. For k > 3 let Wk be the number of non-backtracking closed walks of 
length k in K4, and define = Wk/2k. Further, let Yk be a Poisson random variable 
with expectation Xk, with {Yk}k independent, and define 5k = (— l/2) fc . Then with 
G = K 4 , 



X<! 1 W:=H(l + 6. 



G d . xxr ._ TT f-i , x \Yi g-Aiii 



EX G 



Proof. Let X = X^ 4 . It follows from Examples 13.71 and 14.51 that 



E(X 



2\ 



{Exy 



2 10 3-3/2 5 -l n -3/2_ 



By comparing with Example 15.51 we find that (A4) of [Til Theorem 9.12] is satisfied: 
that is, 

EX 2 \ 
( EX y ~* ex P I 2^ Xk °k J as n ^ 00. 

The other conditions of [TTJ Theorem 9.12] hold, as follows from Lemmas 15.11 and 15.21 
Applying (TTJ Theorem 9.12] completes the proof. □ 



The same argument applies for the multigraph with two vertices and three parallel 
edges, this time using Examples I3.8[ 14.61 and 15.61 leading to the following. 
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Corollary 5.8. Recall that K\ denotes the multigraph with two vertices and three parallel 
edges. For k > 2 let Wk be the number of non-backtracking closed walks of length k, and 
define Xk = Wk/2k. Further, let Y& be a Poisson random variable with expectation Xk, 
with {Yk}k independent, and define 5k = (— l/2) fe . Then with G = K\, 

y. OO 

-A> W := J] (1 + 5 2i ) Y2i e~ x ^. 



EX, 



G 



It is immediate that the limiting distribution W satisfies W > (with probability 
1) in both Corollary 15.71 and 15.81 Hence L n (G) a.a.s. has a perfect matching, for both 
G = K4 and G = K\. This also follows from [12J. 



6 Summation by Laplace's method 

In this section we prove our main approximation tool, Theorem 12.31 which performs a 
summation over lattice points. We will require a little more theory about lattices. The 
following surprising duality was proved by McMullen [14J. (See also [19].) 

Lemma 6.1. Let V be a subspace ofW N and let V 1 - be its orthogonal complement. Let 
£ and C be the lattices V H 1< N and V L fl Z N , and assume that the rank of C equals 
the dimension of V (i.e., that £ spans V). Then C has rank dim(K J -) = N — dim(V) 
and 

det(£ ± ) = det(£). 

For our purposes we need a simple extension. 

Lemma 6.2. Let < m < N . Let x\, . . . , x m be linearly independent vectors in 7, N . Let 
V be the subspace ofM. N spanned by x\, . . . , x m and let V 1 - be its orthogonal complement; 
thus 

V ± = {yeR N : (y,Xi) = Ofori = l,...,m}. 

Let C and C 1 - be the lattices V fl Z N and V 1 - fl Z N , and let C be the lattice spanned 
by x\, . . . , x m (i.e., the set {Y^Li n i x i ni E Z} of integer combinations). Then C 1 - has 
rank N — m and 

det(£ x ) = det(£) = det(£ )/q, 

where q is the order of the finite group C/C . Explicitly, q is the number of solutions 
(t u ...,t m ) in (R/Z) rn (or (Q/Z) m ) of the system 

^2xijti = Q (modi), j = l,...,N, (6.1) 

i 

where X{ = {xij)^ =l for i — 1, . . . , m. 
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Proof. Since rank(£) = m = dim(V), we can apply Lemma 16.11 and conclude that 
rank(£ ± ) = N — m and det(£ x ) = det(£). 

Next, Co C V PI = £; moreover, £ and £ both span V and have thus the same 
rank. Hence Lemma [2721 shows that £/£o is finite and det(£) = det(£o)/o. Note further 
that £Cl/ = {^ : U G R} and thus 

g = |£/£ | = e [0,l) m : G £} . 

i 

Furthermore, 

tjXj G £ -<=>- tiXi G Z w •<=>- Xijti = (mod 1) for j = 1, . . . , J, 

i i i 

and the characterization of q follows. □ 
The proof of Theorem 12.31 involves reduction to a special case, which we prove first. 
Lemma 6.3. Suppose the following: 

(i) £ C R r is a lattice with full rank r. 

(ii) K C W is a compact convex set with non-empty interior K° . 

(iii) : K — > R zs a continuous function with a unique maximum at some interior 
point xq G K°. 

(iv) is twice continuously differentiable in a neighbourhood of xq and the Hessian 
H := D 2 <p(xo) is strictly negative definite. 

(v) ip : K\ — > R «s a continuous function on some neighbourhood Ki C K of xq with 
ip(x ) > 0. 

(vi) For eac/i positive integer n there is a vector £ n G R r . 

(vii) For eac/i positive integer n there is a positive real number b n and a function 
a n : (£ + 4) H — > R st/c/i t/iat, as — > oo, 

a„(£) = 0(6 n e™^ + °^), £ G (£ + 4) n ntf, (6.2) 

and 

a n {£) = b n (^(£/n) + o(l))e^ (£/n) , I G (£ + 4) n n#i, (6.3) 

uniformly for £ in the indicated sets. 
Then, as n — > oo, 



£e(c+e n )rmK 



det(£) det(-H) 1/2 



(6.4) 
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Proof. We begin with a few simplifications. We may obviously assume that b n = 1. 
Furthermore, by subtracting (f)(x ) from 0, and dividing a n (£) by e n ^ x °\ we may suppose 
that <fi(xo) = 0. 

Since xq is an interior maximum point, the gradient D<f)(xo) vanishes, and a Taylor 



expansion at xq shows that, using (iv) , as |a: — xq| — ^ 



4>{x) — |(x — so, D 2 cj)(xo)(x — Xq)) + o(\x — xo\ 2 ) (6.5) 

< — C\\x — Xq\ 2 + o(\x — Xq\ 2 ) 

for some positive constant c\. Consequently, there exists 5 > such that the neighbour- 
hood {x : \x — xq\ < 5} is contained in Ki and 

4>{x) < — c 2 \x — x \ 2 , \x — x \ < 5 (6.6) 

for some positive constant ci- We divide the sum in (16.41) into three parts: 

Si := ^ 5 2 := Yl ' Ss := ■ 

|£/n-icol<™~ 1/3 n- 1 / :i <\e/n-x \<8 \i/n-x \>5 

In the sum 5*2 we use f 16 . 3 j) and (j6.6P ; thus each term is 



a, 



,(£) = 0(e n ^ e/n) ) =0(e- C2nl/:i ). 



Since the number of terms is 0(n r ), we obtain S 2 = o(l). 

Similarly, by compactness, if \x— xq\ > 5, then <p(x) < — C3 for some positive constant 
c 3 . Consequently, for large n, (|6.2p shows that each term in S 3 is 

a n (£) = (e n ^ n ^ n / 2 ) = 0(e~ c * n/2 ). 

Again, the number of terms is 0(n r ) and we obtain S3 = o(l). 

We convert the sum S\ into an integral by picking a unit cell U of the lattice £ and 
defining a n (y) := a n (£) for y G U + £, £ G £ + £ n . Let Q n := {j\n /n _ Xo \ <n -i/z{U + f), and 

let Q n := {2; : nx + v^z G Then 

S\ = det(£) _1 / a n (y) dy = det(£) _1 n r//2 / a„(nxo + v^- 2 ) dz. (6.7) 

Note that Q n is roughly a ball of radius n 2 ^ 3 centered at nxo, and Q n is roughly a ball 
of radius n 1 ^ centered at 0. 

If y G Q n , then \y/n — x \ < n^ 1 ^ + 0(n~ l ). Since the gradient D(p(x ) = 0, (iv) 
implies that for x G Q n /n, 

\D(f)(x)\ = 0(\x - x Q \) = 0(n~ 1/3 ). (6.8) 
If y eU + £ C Q n , then |y/n - £/n\ = 0(1 /n) and ([63]) implies 

n<f){y/n) - n<j){£/n) = O^rT^n' 1 ) = 0(n~ 1/3 ), 
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and thus (16.31) implies, uniformly for y G Q n , 

a n {y) = a n (£) = {if>{y/n) + o(l))e n ^ n \ (6.9) 

For every fixed z G W , this and the Taylor expansion (I6.5P show that, as n — > oo, 
using the continuity of if), 

a n {nx + y/Ez) -> iP{x )e^ z ' D2<l)(xo)z) . 

Moreover, (16. 6p and (16.91) provide a uniform bound, for all z G W, 

\a n (nx + v^z)l§ n (z)| < Cie~ C2 ^ 12 . 

Further, lg (2;) — > 1 for every z. Hence, dominated convergence shows that 

a n {nx + Vnz) dz -> / ^(xo)e^ D20(:Co)2) dz 

= ^(x o )(27r)^ 2 det(- J D 2 0(^o)) _1/2 . 



The result follows from this and (16. 7p . together with the estimates S2 = o(l) and 
S3 = o(l) above. □ 

Proof of Theorem \2.3[ First, replacing K by K — w, a n (£) by a' n (£) := a n (£ + nw), £ n 
by ^ n — ww, an d translating <fi and ■?/>, we reduce to the case w = and thus W = V and 

4e V. 

Choose a lattice basis {zi, . . . , z r } of C. Consider the mapping T : IR r — > V C 
given by (2/1, . . . ,y r ) 1— ► 5^i=i which thus maps Z r onto £. We apply Lemma [6731 to 
£' := Z r , K' :=T-\K), (f>oT, ifjoT, £' n := T~ l {£ n ), and a n {T(k)), k G {£' + £' n ) HnK'. 
The Hessian -D 2 (0 o T)(T _1 a;o) equals (H(zi, Zj))^- V and its negative has determinant, 
by ([23]) and (Q, 

det ^))r, i= i = det(-H\ v ) det(( Zi , 2 i »^- =1 = det(-#| y ) det(£) 2 . (6.10) 

Hence, (I2.7P follows from Lemma [6731 Note that the Hessian D 2 (<f)oT)(T~ 1 xo) is always 
negative semi-definite, because xq is a maximum point. Hence, it is negative definite 
unless its determinant is zero, which is ruled out by (I6.10p and the assumption that 
det(-H\ v ) ^ 0. □ 
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